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Abstract. 2-kaon and 2-pion correlation functions for an expanding thermal- 
ized source are compared. In the Yano-Koonin-Podgoretskii parametrization 
7—i of the correlation function, the HBT radius parameters are shown to obey Mj_- 

£H scaling in the absence of collective transverse flow. This scaling is broken by 

^ transverse flow. An accurate comparison of pion and kaon correlations can 

^"5 thus resolve issue whether the observed M±_ -dependence of the transverse ra- 

dius parameter is due to transverse collective flow or other transverse gradients. 
(N Effects from resonance decays are shortly discussed. 

^T 1 1. Introduction 

In the last few years a large body of evidence has been accumulated that the hot 
and dense collision region in ultrarelativistic heavy ion collisions thermalizes and 
\Q shows collective dynamical behaviour. This evidence is based on a comprehensive 

ON analysis of the hadronic single particle spectra. It was shown that all available 

data on hadron production in heavy ion collisions at the AGS and the SPS can be 
Y understood within a simple model which assumes locally thermalized momentum 

r ^ distributions at freeze-out, superimposed by collective hydrodynamical expansion 

^ in both the longitudinal and transverse directions [l]-[6]. The collective dynamical 

f| behaviour in the transverse direction is reflected by a characteristic dependence of 

the inverse slope parameters of the mj_ spectra ("effective temperatures") at small 
mj_ on the hadron masses [1, 2]. New data from the larger Au+Au and Pb+Pb 
systems [7, 8] support this picture and show that the transverse collective dynamics 
is much more strongly exhibited in larger collision systems than in the smaller ones 
from the first rounds of experiments. The amount of transverse flow appears also 
to increase monotonically with collision energy from GSI/SIS to AGS energies, but 
may show signs of saturation at the even higher SPS energy [9]. 

It must be noted, however, that the "proof of thermal and hydrodynamic be- 
haviour from an analysis of single particle spectra strongly relies on model assump- 
tions. A model independent extraction of flow velocities and thermal freeze-out 
temperatures from the measured momentum spectra is not possible [2, 4]. Ex- 
plicit dynamical model calculations and additional theoretical consistency argu- 
ments [3, 10, 11] are required to achieve this goal, and there have been numerous 
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alternative suggestions to explain the observed features of the hadron spectra with- 
out invoking hydrodynamical flow, but using differently motivated parametrizations 
[12]. The reason for this ambiguity is that single particle momentum spectra, as a 
matter of principle, contain no direct information on the space-time structure and 
the space-momentum correlations induced by collective flow [13]: In terms of the 
phase space density at freeze-out ("emission function") S(x,p) the single-particle 
spectrum is given by EdN/d 3 p = f d 4 x S(x,p), and one sees that the space-time 
information contained in S is completely washed out by integration. Thus, on 
the single-particle level, one must perform comprehensive model analyses to show 
that a simple hydrodynamical model with only a few thermodynamic and collec- 
tive parameters can fit all the data, and one must use additional consistency checks 
in order to show that the values of the extracted fit parameters lead to a sensi- 
ble and internally consistent theoretical picture. The published literature abounds 
with examples demonstrating that without such consistency checks the theoretical 
ambiguity of the single particle spectra is nearly infinite. 

This is the point where two-particle Bose-Einstein correlations between the 
momenta of pairs of identical particles provide crucial new input. They give direct 
access to the space-time structure of the source and its collective dynamics, and 
although some remaining model dependence cannot be avoided the set of possible 
model sources is drastically reduced. The two-particle correlation function C(q, K) 
is usually well approximated by a Gaussian in the relative momentum q, with 
Gaussian width parameters ("HBT (Hanbury-Brown/Twiss) radii") which depend 
on the total momentum A' of the particle pair. It was recently shown [14]-[17] that 
these radius parameters measure certain combinations of the second order variances 
of the space-time structure of the source (the so-called "regions of homogeneity" 
[18] of the source around its point of maximum emission). In general they mix the 
spatial and temporal structure of the source in a nontrivial way [15, 16], and the 
mentioned model dependence enters when trying to unfold these two aspects. 

Collective dynamics of the source leads to a dependence of the HBT radii on 
the pair momentum K; this has been known for many years [19, 20, 21], but was 
recently quantitatively reanalyzed, both analytically [22, 16, 17, 23, 18] and numer- 
ically [24, 13, 25]. The velocity gradients associated with collective expansion lead 
to a dynamical decoupling of different source regions in the correlation function, 
and the HBT radii measure the size of the resulting space-time regions of homo- 
geneity of the source [21, 18] around the point of maximum emissivity for particles 
with the measured momentum K. The velocity gradients are smeared out by a 
thermal smearing factor arising from the random motion of the emitters around 
the fluid velocity [16]. Due to the exponential decrease of the Maxwell-Boltzmann 
distribution, this smearing factor shrinks with increasing transverse momentum K±_ 
of the pair, which is the basic reason for the A'j_-dependence of the HBT radii. 

Unfortunately, other transverse and longitudinal gradients in the source (for 
example spatial and temporal temperature gradients [22, 23]) can also generate a 
A'-dependence of the HBT radii [16]. Furthermore, the pion momentum spectra 
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in particular are affected by resonance decay contributions, but only at small K±_ . 
This may also affect the HBT radii in a A'j_-dependent way [26]. The isolation of 
collective flow, in particular transverse flow, from the K j_ -dependence of the HBT 
radii thus requires a careful study of these different effects. 

We will here study this A'-dependence of the HBT radius parameters within a 
simple analytical model for a finite thermalized source which expands both longitu- 
dinally and transversally. We will use the Yano-Koonin-Podgoretskii parametriza- 
tion [27, 28, 17, 13] of the correlation function which, for sources with dominant 
longitudinal expansion, was shown to provide an optimal separation of the spatial 
and temporal aspects of the source [17, 13]. We point out that, at least in the 
absence of resonance decay contributions, it is a generic feature of all such ther- 
mal models that without transverse collective flow the YKP size parameters show 
perfect Mj_-scaling, i.e. do not depend explicitly on the rest mass of the particles. 
Only the transverse gradients induced by a nonvanishing transverse flow can break 
this Mj_-scaling. Thus accurate measurements of the YKP correlation radii with 
both pions and kaons can identify the presence of transverse collective flow in a 
rather model-independent way. 

At the present moment it is still not absolutely clear to what extent this 
statement must be qualified when resonance decays are included. We have checked, 
however, that the K j_ -independence of the the transverse HBT radius R±_ in the 
absence of transverse flow and other transverse gradients is not modified by res- 
onance decays [29]. Thus, if there is no transverse flow, resonance decays cannot 
fake a transverse flow signal 

2. The source model 

For our study we use the model of Ref. [17] for an expanding thermalized source: 

K-u(x) (t-tq) 2 r 2 (??-T?o) 2 j 
T(x) 2(Ar) 2 2R 2 2(Ar?) 2 J " U 

Here r 2 = x 2 +y 2 , the spacetime rapidity f] = \ ln[(t+z)/(t— z)] and the longitudinal 
proper time r = \/t 2 — z 2 parametrize the spacetime coordinates x 1 * , with measure 
d A x = Tdrd^rdrdcj). Y = ±ln[(E K + K L )j{E K - K L )\ and M L = ^Jm 2 + K\ 
parametrize the longitudinal and transverse components of the pair momentum 
K. T(x) is the freeze-out temperature, R is the transverse geometric (Gaussian) 
radius of the source, To its average freeze-out proper time, At the mean proper 
time duration of particle emission, and Ar] parametrizes [16] the finite longitudinal 
extension of the source. The expansion flow velocity u^(x) is parametrized as 

u^(x) = (cosh rj cosh r/t(r), sinh r/t(r) e r , sinh r\ cosh r/t(r)) , (2) 



(27t)V27t(At) 2 
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with a boost-invariant longitudinal flow rapidity r/i = r\ {v\ = z/t) and a linear 
transverse flow rapidity profile 

m{r) = Vf ("^) • (3) 

r)f scales the strength of the transverse flow. The exponent of the Boltzmann factor 
in (1) can then be written as 

K ■ u(x) = Mj_ cosh(Y — rj) cosh r/t(r) — K±_ -e r sinh r/t(r) . (4) 

For vanishing transverse flow (r/f =0) the source depends only on Mj_. 

3. HBT radius parameters 

From the source function (1) the correlation function is calculated via the relation 

Jd 4 xS(x, K)e^< s -P^ 



2 

\Jd 4 xS(x, A)e 8 '« 



|J d 4 x S(x, K) | 2 |J d 4 x S(x, K) | 2 



(5) 



(see [19, 30-32]). Here q = pi — P2, K = (pi + P'i)l^ are the relative and average 
4-momenta of the boson pair. The quality of the approximation in (5) is discussed 
in [16]. Since pi,P2 are on-shell and thus K ■ q = 0, q° = (3-q (with (3 = K / K° k, 
K/Ek) is not an independent variable and can be eliminated (second equality in 
(5)). Therefore the Fourier transform in (5) cannot be inverted without a model 
for S(x,K). This is the reason for the model dependence of the interpretation of 
HBT correlation data mentioned in the Introduction. 

We use a cartesian coordinate system with the z-axis along the beam direc- 
tion and the «-axis along A'j_ . Then [3 = 0, /?;). We assume an azimuthally 
symmetric event sample and write the C(q, K) in the YKP parametrization [17, 13]: 



C(K,q) = 1 + exp 



-«l9i-4 (9?-(9°) 2 )-(^ + 4) (9-£0 2 



(6) 



Here q\ = qj. + q'y, and i?y, i?o, U are four A'-dependent parameter functions. 
U(K) is a 4-velocity with only a longitudinal spatial component: 

U(K)= 7 (A) (l,0,0, KA")) > Wlth ^ = -/f=f- ( ? ) 

Its value depends, of course, on the measurement frame. The "Yano-Koonin veloc- 
ity" v(K) can be calculated in an arbitrary reference frame from the second order 
space-time variances of the source S(x,K) (for explicit expressions see Ref. [13]). 
It is, to a good approximation, the longitudinal velocity of the fluid element from 
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which most of the particles with momentum K are emitted [17, 13]. For sources 
with boost-invariant longitudinal expansion velocity the YK-rapidity associated 
with v(K) is linearly related to the pair rapidity Y [13]. 

The other three YKP parameters do not depend on the longitudinal velocity of 
the observer system. Their physical interpretation is easiest in terms of coordinates 
measured in the frame where v(K) vanishes. There they are given by [17] 

Rl(K) = if), (8) 
R\(K) = ((z-(A/^)£) 2 )-(A/^j.) 2 (y 2 )«(z 2 ), (9) 
Rl(K) = ^i-i/p L ) 2 )-{f)/pl^{P). (10) 

Here (f(x)) = f d 4 x f(x) S(x, K)j f d 4 x S(x, K) denotes the (A'-dependent) aver- 
age over the source function S(x,K), and x = x — x(K) etc., where x(K) is the 
source point with the highest intensity at momentum A'. Aj_, Ay and Ao thus 
measure, approximately, the (A'-dependent) transverse, longitudinal and temporal 
regions of homogeneity of the source in the local comoving frame of the emitter. 
The approximation in (9,10) consists of dropping terms which vanish in the absence 
of transverse flow and were found in [17] to be generically small even for finite trans- 
verse flow (see below). Note that it leads to a complete separation of the spatial 
and temporal aspects of the source. 

Since in the absence of transverse flow the /3-dependent terms in (9) and (10) 
vanish and the source itself depends only on M±_ , it is clear that then all three 
YKP radius parameters show perfect Mj_-scaling. Plotted as a function of M±_ , 
they coincide for pion and kaon pairs. For nonvanishing transverse flow this Mj_- 
scaling is broken by two effects: (1) The exponent of the thermal factor in (1) 
receives an additional contribution proportional to A'j_ = yj AA? — m 2 . (2) The 
terms which were neglected in the second equalities of (9,10) are non-zero, and 
they also depend on = K_l/Ek- Both effects induce an explicit dependence on 
the particle rest mass and destroy the Mj_-scaling of the YKP size parameters. 

4. Results 

In Fig. 1 we show numerical results for the YKP correlation radii for the source 
(1), with the parameters R = 3 fm, To = 3 fm/c, At = 1 fm/c, Ar] = 1.2, and a 
constant temperature (no temperature gradients!) T = 140 MeV. The dashed lines 
refer to pion, the solid lines to kaon correlations, both at pair rapidity Y = in the 
source cm. frame. Resonance decays are not included (but see discussion above). 

In the left column we show the HBT radii for a source without transverse 
expansion. Due to the absence of transverse temperature gradients, the transverse 
HBT radius Aj_ is independent of M±_ and equal to the geometric Gaussian radius 
R of the source. The strong Mj_-dependence of the longitudinal radius Ay and the 
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effective lifetime Rq is caused by the strong longitudinal expansion with i>/ = z/t. 



pion pairs 
kaon pairs 



■■ 0.6 




m^OO 400 m K 600 800 1000n\200 400m K 600 800 1000 
M ± (MeV/c 2 ) M ± (MeV/c 2 ) 

Fig.l. The YKP radii iJy, and Ro, for pairs with Y = as 
functions of Mj_ = (m 2 + A'^) 1 / 2 . Dashed lines: pion correlations 
(m = m^); solid lines: kaon correlations (m = ttik)- Left column: 
no transverse expansion of source; right column: transverse ex- 
pansion according to (3) with rjj = 0.6. For more discussion see 
text. 

For large Mj_, Ro approaches the width At of the r-distribution in (1), but for small 
Mj_ it is significantly larger. Over the finite longitudinal region i?y from which the 
correlator receives contributions, the proper time hyperbola at r=ro covers a finite 
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range At (in any fixed frame, in particular in the source rest frame). i?o measures 
this At plus the Gaussian width At. Since i?y decreases for large Mj_, so does At, 
and only At survives in the limit. 

For vanishing transverse flow rjj = 0, all three YKP radii show perfect Mj_- 
scaling. If the source expands longitudinally like ours, i?y and i?o are thus smaller 
for kaon pairs than for low-momentum pion pairs, even without resonance decays. 
Due to the smaller i?y for kaons, they cannot probe the longitudinal variation At 
of time along the freeze-out surface, and i?o for kaons essentially measures only the 
Gaussian widths At of the proper emission time distribution. 

The right column corresponds to a transversally expanding source with r/f=0.6. 
This value is probably at the upper edge of the range required to explain present 
heavy ion data [4, 7, 33]. One sees three types of effects: (a) Now also the trans- 
verse radius develops an Mj_-dependence. The M±_ slope is directly related to 
the magnitude of rjj [29, 34]. (b) The Mj_-scaling of the YKP radii is broken. The 
effect is not very big, but it goes in opposite directions for i?y and the two other 
radius parameters. This will be helpful in disentangling flow from possible reso- 
nance decay effect which are stronger for pions than for kaons and always tend to 
increase all three radius parameters, (c) i?o receives contributions from the last two 
terms in the middle expression in Eq. (10) and no longer just measures the effective 
lifetime \/ {t 2 ) . (The corresponding contributions to i?y in Eq. (9) still vanish, due 
to z-symmetry at Y = 0.) The scale-breaking effects from the correction terms are 
in fact stronger and of opposite sign relative to the lifetime \/ (t 2 ) . The relative 
difference between i?o and \/ (t 2 ) is expected to decrease for larger systems with 
larger values of To and At. In this sense i?o is still a good measure for the effective 
lifetime, actually more so for kaons than for pions. 
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